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Abstract 

The proper-time resolution function for the measurement of the time evolution of B 
mesons with the Belle detector at KEKB is studied in detail. The obtained resolution 
function is applied to the measurement of B meson lifetimes, the B°B° oscillation 
frequency and time-dependent CP asymmetries. 
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1 Introduction 



KEKB is an asymmetric electron-positron collider designed to produce boosted 
B mesons [1]. At KEKB, electrons (8.0 GeV) and positrons (3.5 GeV) collide at 
a small (±11 mrad) crossing angle. Their annihilations produce T(45) mesons 
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moving nearly along the z axis, defined as anti-parallel to the positron beam 
direction, with a Lorentz boost factor of (Prfr = 0.425, and decaying to B°B° 
or B + B~ . Precise determination of the proper-time interval (At) between the 
two B meson decays is essential for the measurement of B meson lifetimes, 
the B°B° oscillation frequency, and time- dependent CP asymmetries. Since 
the B mesons are nearly at rest in the T(4S) center of mass system (cms), At 
can be determined from the separation in z between the two B decay vertices 
(Az). The average Az at KEKB is cr B (/37) T ~ 200 fim, where r B is the B 
meson lifetime. In this article, we consider the analysis in which one of the 
B decay vertices is determined from a fully reconstructed B meson, while the 
other is determined using the rest of the tracks in the event. 

In order to extract the true At distribution from the observed Az distribution 
it is necessary to unfold the vertex detector resolution and (possible) bias in the 
measurement of Az. Because the detector resolution is of the same order as the 
average Az at KEKB, an understanding of the resolution is crucial for precise 
measurements. For measurement of time-dependent quantities, we employ an 
unbinned maximum likelihood method. A probability density function for the 
likelihood fit is obtained as a convolution of a theoretical At distribution with 
the resolution function. 

This paper describes the resolution function developed for the precise measure- 
ment of B meson lifetimes with fully-reconstructed hadronic decay final states 
at the Belle experiment [2]. The obtained resolution function is also applied to 
the measurements of the B°B° oscillation frequency [3] and time-dependent 
CP asymmetries [4]. 

The organization of the paper is as follows: We give a brief description of the 
Belle detector in Section 2. In Section 3, the method of vertex reconstruction 
is described. Details of the resolution function and its application to the data 
are described in Sections 4 and 5, respectively, followed by a conclusion in 
Section 6. 



2 The detector 

The Belle detector, shown schematically in Fig. 1, is a large-solid- angle mag- 
netic spectrometer that consists of a three-layer silicon vertex detector (SVD), 
a 50-layer central drift chamber, a mosaic of aerogel threshold Cherenkov 
counters, time-of-flight scintillation counters, and an array of CsI(Tl) crystals 
located inside a superconducting solenoid that provides a 1.5 T magnetic field. 
An iron flux return located outside of the coil is instrumented to detect K\ 
mesons and to identify muons. The detector is described in detail elsewhere [5]. 
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Fig. 1. Belle detector. 



The SVD, which plays an essential role in reconstructing decay vertices, con- 
sists of three concentric cylindrical layers of double-sided silicon strip detec- 
tors. It covers the angular range of 23° < 9 < 139°, where 6 is the polar angle 
from the z axis. This coverage corresponds to 86% of the full solid angle in 
the cms. The three layers are located at the radii of 3.0, 4.5, and 6.0 cm. The 
strip pitches are 84 /im for the measurement of the z coordinate and 25 fim for 
the measurement of azimuthal angle <fi. The impact parameter resolutions for 

charged tracks [6] are measured to be a xy = J (19) 2 + (50/(p/3 sin 3//2 6)) 2 fim in 



the plane perpendicular to the z axis and a z = V(36) 2 + (42/(p/3sin 5 / 2 9)) 2 
along the z axis, where f3 = pc/E, p and E are the momentum (GeV/c) and 
energy (GeV) of the particle. 



3 Proper-time interval reconstruction 

In this section, we describe the reconstruction of the proper-time interval 
between the decay points of the two B mesons produced at KEKB. Figure 
2 illustrates reconstruction of the decay vertices. The decay vertices of the 
two B mesons in each event are fitted using tracks that have at least one 
three-dimensional coordinate determined from associated r-<f> and z hits in 
the same SVD layer plus one or more additional z hits in other SVD layers. 
We impose the constraint that they are consistent with the interaction point 
(IP) profile, smeared in the r-(f> plane by 21 /im to account for the transverse B 
decay length. The IP profile is described as a three-dimensional Gaussian, the 
parameters of which are determined in each run (in finer subdivisions, every 
10,000 - 60, 000 events, for the mean position) using hadronic events. The size 
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of the IP region is typically a x ~ 100 //m, a y ~ 5 /im, and cr 2 ~ 3 mm, where 
x and |/ denote the horizontal and vertical directions, respectively. 

One B meson is fully reconstructed in one of the following decay modes 2 : 
B° -> D+n-, D* + ti-, D*+p~, J/^K° S , J/tpK* , B~ -> D°tt~, and J/^K~, 
where J/^ is reconstructed via J/^ — > £ + £~(£ = e,fi) decay, and -D* + via 
D* + — > Z}°7T + decay. Neutral and charged D mesons are reconstructed in the 
following channels: D° — > K~ir + , fT~7r + 7r°, f^~7r + 7r + 7r~, and _D + — > K~ir + ir + . 
The details of the event selection can be found elsewhere [2]. 

In the case of a fully reconstructed B — > D^X decay, the B decay point is 
obtained from the vertex position and momentum vector of the reconstructed 
D meson and tracks other than the slow n + candidate from D* + decay. For a 
fully reconstructed B — > J/ipX decay, the B vertex is determined using lepton 
tracks from the J/ip. 

The decay vertex of the associated B meson is determined by applying a 
vertex-fit program to all tracks not assigned to the fully reconstructed B me- 
son; however, poorly reconstructed tracks (with a longitudinal position error 
in excess of 500 pm) as well as tracks that are likely to come from K s decays 
(forming the Kg mass with another track, or emanating from a point more 
than 500 /im away from the fully reconstructed B vertex in the r-<f> plane) are 
not used. If the reduced x 2 (xV n -d.f.) associated with a found vertex exceeds 
20, the track making the largest x 2 contribution is removed and the vertex is 
refitted. This procedure is repeated until x 2 /n.d.f. < 20 is obtained or only 
one track is left. If, however, the track to be removed is a lepton with a cms 
momentum greater than 1.1 GeV/c, we keep the lepton and remove the track 
making the second largest x 2 contribution. This is because high-momentum 
leptons are likely to come from primary semi-leptonic B decays. The presence 
of a secondary charm (b — > c) decay vertex in the associated B meson results 
in a shift of the reconstructed vertex point toward the charm flight direction 
and degrades the vertex resolution. A Monte Carlo (MC) simulation study 
shows that the shift and the resolution of the associated B decay vertex are 
~ 20 /zm and ~ 140 fj,m (rms), respectively, while the resolution of the fully 
reconstructed B decay vertex is ~ 75 /im (rms). 

Because of the IP profile constraint it is possible to reconstruct a decay vertex 
even with a single track. The fraction of single-track vertices is ~ 10% and 
~ 22% for the fully reconstructed and associated B decays, respectively. For 
the multiple-track vertices, the quality of the vertex fit, for both the fully 
reconstructed and associated B meson decays, is further evaluated. Using MC 
simulation, we find that the vertex-fit x 2 is correlated with the B decay length 
due to the tight IP constraint in the transverse plane. To avoid this correlation, 

2 Throughout this paper, when a decay mode is quoted the inclusion of charge 
conjugate mode is implied. 
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Fig. 2. Illustration of vertex reconstruction of two B decay vertices. 

we use the variable based on the z information only (contrary to the vertex-fit 
X 2 calculated using three dimensional information): 



n 

£ = (1/2™) ^2 [(Xftcr ~~ Z before)/ £ before 



where n is the number of tracks used in the fit 3 , z before and z* fter are the z 
positions of each track (at the closest approach to the origin) before and after 
the vertex fit, respectively, and £ before is the error of -2 before . Figure 3 shows 
the £ distributions for the (a) fully reconstructed and (b) associated B decay 
vertices, obtained using a MC simulation. Because the x 2 / n -d.f . requirement is 
imposed on the associated B decay vertices, the £ distribution for those vertices 
does not show an extended tail compared to that for the fully reconstructed B 
decay vertices. Figure 4 shows the £ distributions as a function of the B decay 
length for the (a) fully reconstructed and (b) associated B decay vertices. As 
indicated, £ does not depend on the B decay length. We require £ < 100 for 
both vertices to eliminate poorly reconstructed vertices. We find that about 
3% of the fully reconstructed and 1% of the associated B decay vertices are 
rejected in the data. 

The proper-time interval between the fully-reconstructed and the associated 
B decays is calculated as 



At = (z h 



ill 



(2) 



where z lVL \ and z asc are the z coordinates of the fully-reconstructed and associ- 
ated B decay vertices, respectively. We reject a small fraction (~0.2%) of the 
events by requiring |Ai| < 70 ps (~45tb). 



For single-track vertices £ cannot be defined. 
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Fig. 3. The £ distributions for the (a) fully reconstructed and (b) associated B 
decays, obtained using B° — ► J/t/jK® MC events. 
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Fig. 4. Monte Carlo £ distribution as a function of the B decay length, obtained for 
the (a) fully reconstructed and (b) associated B decay vertices. 



4 Resolution function 



4-1 Overview 



We extract the lifetimes of B mesons using an unbinned maximum likelihood 
fit to the observed At distributions. We maximize the likelihood function 
L = YiiP(Ati), where P(Atj) is the probability density function (PDF) for 
the At of the i-th event and the product is taken over all the selected events. 

The function P(At), expressed as 

P(At) = (1 - /oi) [isi g Psi g (At) + (1 - / sig )P bkg (At)] + /oiPoi(At), (3) 

contains contributions from the signal and the background (P S i g and Pbkg), 
where / S i g is the signal purity determined on an event-by-event basis, P S i g is 
described as the convolution of a true PDF (P S i g ) with a resolution function 
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(i?si g ) and Pbkg is expressed in a similar way: 

/+oo 
d(At , )Psig(bk g ) {At')R sig(hkg) (At - At'). (4) 
-oo 

To account for a small number of events that give large At in both the signal 
and background (outlier components), we introduce a fraction of outliers (/ Q i) 
and a Gaussian function (P oi (At)) to model its distribution. The true PDF 
for the signal, P sig (At; t#), is given by 

V sig (At; t b ) = ^~ exp ( -1M ] , (5) 

where r B is, depending on the reconstructed mode in the event, either the 
B° or the B~ lifetime. The signal PDFs for the measurements of the B°B° 
oscillation frequency and the time-dependent CP asymmetry are given in 
Section 5. 

The resolution function of the signal is constructed as the convolution of four 
different contributions: the detector resolution on Zf u \ and z asc (Rf n \ and -Rase); 
an additional smearing on ^ asc due to the inclusion of tracks which do not 
originate from the associated B vertex (i? np ), mostly due to charm and Kg 
decays, and the kinematic approximation that the B mesons are at rest in the 
cms (i?k)- The overall resolution function, R sig (At), is expressed as 

R sig (At) = f[[ +0 ° d(At') d(Af') d(Af) R ful (At - At^R^At' - At") 

J J J — 00 

xR np (At" - At"')R k (At" f ). (6) 

We use a MC simulation to understand the resolution function and determine 
its functional form. The MC events are generated using the QQ event genera- 
tor [7] and the response of the Belle detector is modeled by a GEANT3-based 
full-simulation program [8]. One of B mesons in each event is forced to de- 
cay to J/il)K° S) J/ipK' or Dir while the other decays generically to one of all 
possible final states. 



4-2 Detector resolution 

The detector resolution (-Rf u i and -Rase) is studied using a special MC simula- 
tion in which all short-lived (r < 10~ 9 s) secondary particles (including K s 
and A) are forced to decay with zero lifetime at the B meson decay points. 
Figures 5 (a) and (b) show the distributions of the difference in z between the 
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(a) fully reconstructed B 



(b) associated B 
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Fig. 5. The 5z distributions of the (a) fully reconstructed and (b) associated B 
vertices obtained from a B° — > J '/i^l 
of a fit to the sum of two Gaussians. 

reconstructed and generated vertex positions: 



5z a 



gen 



(7) 



where q = ful(asc) is for the fully reconstructed (associated) B vertex, and 
the superscripts 'rec' and 'gen' denote the reconstructed and generated vertex 
positions, respectively. Results of the fit to a sum of two Gaussians are also 
shown. The fitted curves do not represent the 5z distributions in the tail 
regions. We also find that even a sum of three or more Gaussians with constant 
standard deviations cannot represent 5z properly. We therefore consider a 
more elaborate function that uses the vertex-by-vertex z-coordinate error of 
the reconstructed vertex, o~ z q (q = full, asc), as an input parameter. The value 
of a z q is computed from the error matrix of the tracks used in the vertex fit 
and the size of the IP region. To construct functional forms of i?f u i and i? asc 
we investigate the distribution of a pull, defined as Sz q divided by a z q . If the 
<7g estimation is correct on average, the pull distribution is expected to be a 
single Gaussian with the standard deviation of unity. 

Because the resolution for the multiple-track vertices is better than that for 
the single-track vertices, we consider them separately. Figure 6 shows the dis- 
tributions of of ull and a z sc for the multi-track and single-track vertices obtained 
from B° -> Jji)K% data. 



4-2.1 Multiple-track vertex 

We investigate the vertex fit quality dependence of the resolution using the 
value of £ defined in Eq. (1). We find that a pull distribution for vertices with 
similar £ values can be expressed as a single Gaussian. Furthermore, we find 
that the standard deviation of the distribution has a linear dependence on £ 
as shown in Fig. 7. 
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Fig. 6. Distributions of the vertex-by-vertex z-coordinate error a z of the (a) (b) 
mutli-track and (c) (d) single-track vertices for the fully reconstructed and associ- 
ated B decays, obtained from B° — ► J/ipKg data. 




Fig. 7. Standard deviations of the pull distributions as a function of £ for (a) the 
fully reconstructed and (b) the associated B meson vertices. The distributions are 
obtained from aB -* J/iJjK^ MC sample. 
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Fig. 8. Distributions of (a) 5zf u \ and (b) 5z a sc for multiple-track vertices, with 
i?^ tiple (fefui) and -RSc ltiplc (<J-Zasc)) respectively. Figure (c) is the S(Az) distribution, 
and the convolution of R^ lUple {5z{ u i) and R^ ltiple (Sz asc ). 

Results from this MC study lead us to model the detector resolution of the 
multiple-track vertex using the following function: 

R^(5z q ) = G(5z q , (s° q + s\i)al) (q = ful, asc), (8) 

where G is the Gaussian function, 



G(x] a) 



27nx 



exp 



x 
'2^ 



(9) 



The scale factors s q and s q are treated as free parameters and determined 
from the lifetime fit to the data. Figures 8 (a) and (b) show the <5zf u i and 5z asc 
distributions, respectively. Superimposed are the curves obtained by summing 
vertex-by-vertex R q nnltiple (5z q ) functions. The curves well reproduce the 5z q 
distributions. This demonstrates that /£™ ultl P le represents the detector resolu- 
tion better than the sum of two Gaussians in Fig. 5. Figure 8 (c) shows the 
distribution of the residual of Az, 5(Az) = Az vcc — Az gcn together with the 
convolution of /C ltiplc (5 z M ) and #™ ltiplc (<^ asc ). 



4-2.2 Single-track vertex 

For the single-track vertices, £ is not available. The resolution function of the 
single-track vertices, R s q msle (Sz q ) (q = ful, asc), is expressed as a sum of two 
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(a) fully reconstructed B 



(b) associated B 
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Fig. 9. The (a) 5zf u \ and (b) fcasc distributions for single-track vertices with R s ^ g ° 
and i?asc gle . 

Gaussians, one for the main part of the detector resolution and the other for 
the tail part, which is due to poorly reconstructed tracks: 

Rf sle (5z q ) = (1 - f tail )G(5z q ; s main cr*) + f tail G(5z q , s tail a* q ), (10) 

where s ma i n and Staii are global scale factors which are common to all single- 
track vertices. Figure 9 shows the residual distributions of the single-track Zf u \ 
and z asc vertices, together with a fit to R s q mgle (5z q ). 



4-3 Smearing due to non-primary tracks 



We introduce another resolution function, R np , to represent the smearing of 
z^ due to tracks that do not originate from the associated B vertex consists of 
a prompt component, expressed by Dirac's 5-function 6 Daac {8z asc ), and compo- 
nents that account for smearing due to Kg and charm decays. The functional 
form of the non-prompt components is determined from the difference between 
<£asc obtained for the nominal MC sample and that for the special MC sample 
in which all short-lived secondary particles are forced to decay with zero life- 
time at the B decay points, shown in Fig. 10. It can be expressed by a function 
defined as /p-E p (5z asc ; r^ p ) + (1 - f^E^Sz^; r° p ), where / p is a fraction of the 
#2asc > component and E p and E n are: 



E p (x\ t) = — exp ( — ^ for x > 0, otherwise 0, (11) 
E n (x; r) = — exp ^+~^) f° r x — 0> otherwise 0. (12) 
Thus, _R np is given by 

i?n P (^ asc ) = /^ Dirac (^asc) + (l-/« 5 ) [f P E p (5z asc ; r n p p ) + (1 - f p )E a (5 Zasc ; r n n p 

(13) 
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Fig. 11. Average shift of vertex position (z a 



z. 



noNP^ 



noNP 



vs (a) crf sc and (b) £ 

a 



is obtained from a MC sample in which secondary decays are turned off. The 
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events with z a 



< 
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where fs is the prompt- component fraction. We find that the vertex position 
shift has linear dependence on both <rf sc and £f sc as shown in Fig. 11 for multi- 
track vertices. Consequently, we assume r p p and r^ D are bilinear with cr asc and 

£asc &S 



np 



P " -° + rl(sl c + 4c£ a 2 scKsc/c(/37)r, and (14) 



Til = T 



np ' p 1 ' p 



C = r n ° + r^l + sLUeXJcmr- (15) 



We determine six parameters in R np , fs, / p , r° r°, and by fitting 
the convolution of i?™" ltiplc and R np to the 5z asc distributions for B° and _B~ 
separately, as shown in Fig. 12. In this fit, the scale parameters, s° sc and s^ c , 
for i?™" ltiple are fixed to the values (common to B° and 5~) obtained by fitting 
i?™" ltiple to the special MC sample (in which all short-lived secondary particles 
are forced to decay promptly at the B meson decay points). Results, shown 
superimposed, well represent the distributions. 
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Fig. 12. The fcasc distributions of multiple-track vertices for (a) B° — » J/t/jKg and 
(b) B~ -> J/^if- decays. 
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Fig. 13. Average shift of vertex positions (z asc — z 
vertices. The events where z asr < z^?^ p are excluded. 



versus cr asc for single-track 



For single-track vertices we consider the correlation between the vertex posi- 
tion shift and a* sc . Figure 13 shows the vertex position shift versus <rf sc for 
the single-track vertices. 

Since i? asc for the single-track vertices is defined as a sum of main and tail 
Gaussians, we also introduce -R^ am and R^ 1 for main and tail parts, respec- 
tively. Each of -R^ am and R^ 1 is expressed by the function of Eq. (13) with 
parameters defined as: 



( T np)main Tp "F T p s mam cr aS c/ C (/^7) Y 

(Cp)main = T° + T^S^CT^J c(^) T , (16) 

(retail =r p ° + r p 1 s tai i< c /c(/37) T , and 

(C) ta n = r n ° + r n 1 Stail < c /c(/5 7 )T. (17) 



The convolution of R asc and R np for single-track vertices is, thus, defined as: 
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Fig. 14. The 5z asc distributions of single-track vertices for (a) B° and (b) B 
mesons. 



Table 1 

Shape parameters of R np used for the lifetime fit. The values are determined for 
multiple- and single-track vertices, separately, using Monte Carlo fc asc distributions. 
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(18) 



Figure 14 shows the <5z a sc distributions for the single track vertices. The su- 
perimposed curves are the results of a fit to the function given by Eq. (18). 

Table 1 lists the shape parameters of R np determined by fitting -Rasc'Si-Rnp^asc) 
to MC 5z asc distributions. These parameter values are held fixed when the life- 
time fit to the data is performed. 
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4-4 Kinematic approximation 



The proper time interval calculated as Eq. (2) , At = (zf u i — -s; asc )/[c(/37)x], is 
equal to the true proper-time interval when the cms motion of the B mesons 
is neglected. The difference between At and the true proper-time interval 
A^true = tf u i — t asc is calculated from the kinematics of the T(4S') two-body 
decay: 



X = At - Attruc = (^ful - 2 asc )/[c(/?7) T ] - (tfui - t asc ) 

= [tfuic(/?7)ful - 4scC(/37) asc ]/[c(/?7) T ] - (tful - ^asc) 

= mU/iPlh - l]ifui - [(/3tW(/37)t - l]4sc, (19) 

where (Pj)m and (/?7) asc are Lorentz boost factors of the fully reconstructed 
and associated B mesons, respectively, and their ratios to (/?7)r are given as: 



perns cms p n o flcms 

(/57)fui/()97)r = — + ^5 ^~l + 0.165cos#r, and (20) 



pems cms zicms 

()97)«c/(/?7)t = ^- - PB ™" B ~ 1 - 0.165 cos 6^ , (21) 

where /3 T = 0.391 is the velocity of the T(45) in units of c, £ B ms ~ 5.292 GeV, 
p c ™ s ~ 0.340 GeV/c and 6 c B ms are the energy, momentum and polar angle of 
the fully reconstructed B in the cms, and tub is either the B° or the B~ mass. 
The difference x can, therefore, be approximated as 

0.165 cos ^ ms (t ful + t asc ). (22) 

i?k, which accounts for x, can be given as a function of cos# B ms . Because tf m 
and t asc distributions follow E p (t fvl ] r B ) = ^ exp(-t fm /r B ) and £ p (t asc ; r B ) = 
^ exp(— t asc / re), respectively, the probability density of obtaining x and At true 
simultaneously is given by: 



POO P OO 

F(x, At truc ) = J J dt ful dt asc £ p (t ful ; r B )£ p (t asc ; r B )5 Dirac (At truc - (t ful - t 

Xd DiraC (x - {[(/3 7 )ful/(/37)T " l]tful " [(/97)asc/(/97)T - l]*asc 

and the probability density of obtaining At tr uc is given by: 

roo roc 

F ( Attruc ) = / / dtf m dt asc E p (tf m ; T B )E p (t asc ; t b )5 irac (At true — (t fm — t asc )). 
Jo Jo 

(24) 
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Fig. 15. The x = At — At tr ue distribution for B° — ► J/ipKg events together with 
the function R^x). 

Rk(x) is, then, denned as the conditional probability density of obtaining x 
given At tme . It is expressed as R\{x) = F(x, At true )/F(At true ) which gives: 
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' B 

- rn B 



fDirac 



(x-(- 
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f3rm B 
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{cose^ 8 
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Figure 15 shows the x distribution for B° — > J/i[)Kg events with the function 
Rk(x). The expected theoretical At distribution -P(At) can be expressed as a 
convolution of the true PDF P S i g (At true ; r B ) with R k (At — At true ): 

+ for At > 
- for At < 

(26) 



P(At) 



mB ( | At | 

\ y ms prm B 1 " 



4-5 Background distribution 



The signal purity / S i g in Eq. (3) is determined on an event-by-event basis as a 
function of two kinematic variables, the energy difference AE = E B ms — E^ m 
and the beam-energy constrained mass Mt, c = \J (E^ m ) 2 — (p§" s ) 2 , where 
^beam i s the beam energy in the cms. Typical / sig values for the decay modes 
used for the measurement of B meson lifetimes and B°B° oscillation frequency 
are listed in Table 2. The values listed are obtained for the signal region defined 
as |Mb c — m B \ < 3cr. The composition of the background is studied by a MC 
simulation and is listed in Table 3. The largest contribution is from cc events 
produced in the continuum, while other continuum events (u, d, s events) and 
combinatorial backgrounds from BB events also contribute. 
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Table 2 

Typical values of the signal purity / s i g for the decay modes used for the measurement 
of B meson lifetimes and B°B° oscillation frequency. 



Decay mode -Bf u i 


Signal purity / sig 




0.861 


B° -► D* + tt- 


0.812 


B° -► D*+p- 


0.699 


B° -► 

5- -► D°TT- 


0.940 
0.948 
0.699 


B- -► J/V'K- 


0.955 



Table 3 

The composition of the background in the decay modes used for the measurement of 
B meson lifetimes and B°B° oscillation frequency. The sources of the backgrounds 
are the continuum production of uu, dd, ss pairs (qq(u, ,d,s)) and cc pairs (qq(c)), 
and combinatorial backgrounds from B + B~ and B°B° events. B — ► J/tftX modes, 
in which the background fraction is very small, are not listed. 



Decay mode Bf u \ 


qq(u,d,s) 


99(c) 


B+B- 


B°B° 


W -► D+it- 


0.25 


0.40 


0.21 


0.14 


B° L>*+vr- 


0.10 


0.48 


0.24 


0.18 


B° -► L>*+p- 


0.10 


0.43 


0.26 


0.21 


B- _> 


0.29 


0.45 


0.17 


0.09 



The background PDF, Pbkg(At), is modeled as a sum of exponential and 
prompt components (Pbkg(At)) convolved with i?bk g (At): 

/+00 
d(At')V hkg (At')R hkg (At - At'), (27) 
-00 

where 

V hkg (At) = f^5»™(At - /if g ) + (1 - f^)—- exp - |At ^ 1 (28) 

^Tbkg \ Tbkg / 

with jj^ ks and /x^ kg being offsets of the distribution, and 

^bkg(At) = (1 - /S)G(Atj iV^ful + ^sc) + AaifG(At; ^if^ful + ^sc)- 

(29) 

Different values are used for s^ in , s}*f , and f^f depending on whether both 
vertices are reconstructed with multiple tracks or not. All parameters in P bkg (At) 
are determined by the fit to the At distribution of the background-enhanced 



17 



control sample (i.e. events in the sideband region of the AE or Mb c distribu- 
tion). Table 4 lists the values obtained separately for each decay mode used 
in the lifetime fit. 



4-6 Outliers 



We find that there still exists a very long tail that cannot be described by 
the resolution functions discussed above. The outlier term is introduced to 
describe this long tail and is represented by a single Gaussian with zero mean 
and event-independent width: 

P ol (At)=G(At,a ol ). (30) 

The global fraction of outliers / Q i (in Eq. (3)) and the a Q \ are left as free 
parameters in the lifetime fit. Different values are obtained for / Q i depend- 
ing on whether both vertices are reconstructed with multiple tracks or not 

(jmultiple^ ^single^ Wg find ^ ^multiple ig ^ ^ 1Q - 3 &nd jingle _ 1Q - 2 



5 Applications 

5.1 B meson lifetimes and B°B° oscillation 



The lifetimes of the B° and B~ mesons are extracted from a 29 fb _1 data 
sample, which contains 31.3 million BB pairs [2]. In the final fit to the events 
in the signal region, we determine simultaneously twelve parameters: the B° 
and B~ lifetimes (r-g , Tb-), four scale factors (s° ul , s} uh s asc , s asc ) for i?^ ltiple 
and i?™" ltiple , three parameters (s ma j n , s ta ii and / tail ) for i^i and R^ c , and 
three parameters (<r ol , f™ ltiple , /oi mgle ) f° r outlier component. Table 5 lists the 
result of the fit. The fit yields: 



t w = 1.554 ± 0.030(stat) ± 0.019(syst) ps, 
t b - = 1-695 ± 0.026(stat) ± 0.015(syst) ps, and 
t b -/t w = 1.091 ± 0.023(stat) ± 0.014(syst). 

The resulting At resolution for the signal is ~ 1.56 ps (rms) for this data sam- 
ple. Figure 16 shows the distributions of At for B° and B~ events in the signal 
region with the fitted curves superimposed. The systematic errors arising from 
the resolution functions are estimated by comparing the results obtained using 
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Table 4 

The background shape parameters obtained by fitting Pbkg(Ai) to the background- 
enhanced control sample. 



(a) B — > J/ijjK. modes 


Parameter 


^0 JlibK°a 


b° -► j/if>K*° b~ -► J/i>K- 


rmamMultiplc 


0.40 ± 0.12 


i no 4- n 9^ 


n 7q -i- n 9^ 

U. I a :x u.^o 


V^tail .'multiple 


9 

j.'±u_2 57 


R Q7+6-39 
D - y ' -2.30 


1 nn+0-64 
1.90_ 037 


C f bkg, > , ■ , 

V J tail Multiple 


n on x n 1 9 

U.Zc/ X U.1Z 


0.03 ± 0.04 


66+ 21 

U.OD_q 2 g 


( fr bkg ) m „l t ^l 
Vj,5 /multiple 


39 ± 18 


0.08 ± 0.08 


0.85 ±0.05 


V* main jingle 


o.96i8:3 


0.82 ±0.15 


1.03 ±0.10 


v*tail jingle 


4.901™ 


O qq + 2.96 

0.00_2.13 


11 2 +4 - 7 


(/taif) single 


1fi+ - 16 

u - io -0.07 


0.09 ± 0.04 


0.05 ±0.03 


(7,5 § ) single 


0-38±8;§| 


0.18 ±0.13 


0.65 ±0.10 


Tbkg (PS) 


0.39 ±0.25 


1.43 ±0.16 


2 14+ - 41 

z - i4l -0.33 


H? S (P^ 


-0.56 ±0.09 


-0.70 ±0.25 


i -0.00 ±0.05 


^r kS (PS) 


0.23 ±0.18 


-0.00 ±0.14 -0.2ll°;p 


(b) i? — > Z?7r modes 


Parameter 


B° -► D+vr- S° D*+vr- S° -► B- -► D°7T- 


( S ma1n)multiplc 


1.03 ±0.04 




0.90 ±0.07 1.02 ±0.02 


( S taif) multiple 


o nq+0.68 


9 qq + 0.39 

z -°°-o.32 


5.191^58 6.27 ±0.36 


(/taif ) multiple 


0.14 ±0.06 


0.67+_° ? 7 


0.13 ±0.04 0.060 ±0.008 


(/<5 ) multiple 


0.39 ± 0.09 


0.38 ± 0.08 


0.22 ±0.07 0.43 ±0.03 


fs bkg 1 ■ , 

V^mainMnglc 


0.73 ±0.07 


0.87 ±0.10 


0.93 ± 0.08 0.79 ± 0.03 


( S taif) single 




4 cu+4.05 
^■ 04 -1.34 


3.52l°; 4 £ 5.99 ± 0.54 


(/taif ) single 


0.12 ±0.04 


0.08 ± 0.05 


0.17 ±0.05 0.09 ±0.01 


(fs § ) single 


0.29 ±0.14 


0.32 ±0.18 


0.19 ±0.10 0.30 ±0.05 


Tbkg (PS) 


1.10 ±0.14 


1.68^ 


0.87 ±0.11 0.98 ±0.05 


S (P^ 


-0.03 ±0.03 


0.00 ± 0.03 


0.11 ±0.07 -0.02 ±0.01 


^r kS (PS) 


0.00 ± 0.08 


-0.11 ±0.14 - 


-0.13 ±0.07 -0.11 ±0.02 



19 



Table 5 

The parameters determined by the lifetime fit. 



Parameters 


Values 


T B° (P S ) 


1.554 ± 0.030 


t b - (ps) 


1.695 ± 0.026 


s ful 


0.809 ±0.148 


s ful 


0.154 ±0.013 


S° 

•'asc 


0.753 ± 0.065 


s 1 

°asc 


0.064 ± 0.005 


^rnain 

«tail 

/tail 

Ool (PS) 


647+ - 074 

u - 04 ' -0.083 

o nn+2.23 
,3 - uu -0.99 

0.0831^1 
36.2±g;g 



^multiple (5.83+|°^ x 1Q -4 



^.single Q 03Q6 ± Q 0Q36 

a different functional form and by varying the resolution parameters within 
their errors. They amount to 0.010 ps, 0.009 ps and 0.006 for r^ , r B - and 
t b-/ t ~b~o, respectively. MC simulation studies show no bias arising from the 
resolution function in the obtained results. 

The same data sample is used to determine the B°B° oscillation frequency 
Arrid from the time evolution of opposite-flavor (OF; B°B°) and same-flavor 
(SF; B°B°, WW) neutral B decays [3]. The signal PDF is defined as the 
convolution of the true PDF, 



V OF (At) = — !— exp ( -1^1 ] [1 + (1 - 2w) cos(Am d At)] and 
4r-g \ r-go J 

V SF {At) = — *— exp [-1^1] [1 - (1 - 2w) cos(Am d At)}, (31) 
4r "s° V t b» J 

with the resolution function i? sig (At). Here w is the probability for an incor- 
rect flavor assignment and determined simultaneously with Am^ by the fit. 
Figure 17 shows the At distributions for OF and SF events with fitted curves 
superimposed. The fit yields 

Am d = 0.528 ± 0.017(stat) ± O.Oll(syst) ps -1 . 

The error arising from the resolution function accounts for 0.009 ps -1 of the 
total systematic error quoted above. 
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At [ps] 

Fig. 16. The At distributions of (a) neutral and (b) charged B meson pairs, with 
fitted curves. The dotted lines represent the sum of the background and outlier 
component, and the dashed lines represent the outlier component. 



5.2 Time- dependent CP asymmetry 



The resolution function obtained above is also used for the measurement of 
mixing- induced CP violation in the neutral B meson system. The Standard 
Model predicts a CP-violating asymmetry in the time-dependent rates for 
B° and B° decays to a common CP eigenstate fcp, where the transition is 
dominated by the b — ► ccs process: 



T(B° - f OP ) - T(B° - /, 



CP) 



T(B»^fcp) + T(BV^f C p) 



-£/ sin 20i sin(Amrft), (32) 



where T(B°,B° — > fcp) is the decay rate for B° or B° to f C p at a proper 
time t after production, £/ is the CP eigenvalue of fcp, and sin 20! is the 
CP violation parameter. A non-zero value for sin20x establishes that CP 
symmetry is violated in the neutral B meson system. We use events in which 
one of the B mesons decays to fcp at time tcp, and the other decays to a 
self-tagging state / tag , which distinguishes B° from B°, at time i tag . The CP 
violation manifests itself as an asymmetry A(At), where At is the proper time 
interval between the two decays: At = t C p — t tag . 



21 




Fig. 17. The At distributions for (a) SF and (b) OF events with the fitted curves 
superimposed. The dashed, dotted, and solid curves show the background, outliers, 
and the sum of backgrounds and signal, respectively. 

The PDF expected for the signal distribution is 

V sig (At,q,w,£ f ) = —^—exp ( — — J [1 - q£ f (l - 2w) sin20i sin(Am d At)], 
4r S o \ Tgo J 

(33) 

where q has the discrete value +1(— 1) when the tag-side B meson is likely to 
be a B (B°), and w is, as mentioned in the previous section, the probability 
for an incorrect flavor assignment (wrong-tag probability) [9]. The PDF is 
convolved with R sig (At) to determine the likelihood value for each event as a 
function of sin20i: 



Pi = {l- /oi)[/si g / V sig (At', q, w, Z f )R*g(At - At')dAt' 

+ (1 - / sig )P bkg (At)] + / iP i(Ai). (34) 

The only free parameter in the final fit is sin20i. The quality of the vertex fit 
for the tag-side B meson (or associated B meson) can be correlated with w, 
and, therefore, the resolution function can also be correlated with w primarily 
because of smearing due to non-primary tracks. R np is designed to account 
for this correlation by making the parameters of vertex position shifts (r^ p 
and r° p) a function of vertex-fit qualities (<r acs and £ aS c)- MC simulation study 
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At (ps) 

Fig. 18. The At distributions for the events with q£f = +1 (solid points) and 
q^f = — 1 (open points). The results of the global fit with sin20i = 0.719 are shown 
as solid and dashed curves, respectively. See text for the detail. 

shows that the remaining w dependence is negligible. 

The value of sin 20! is measured using a 78 fb _1 data sample, which con- 
tains 85 million BB pairs [4]. This data sample has been analyzed using a 
new track reconstruction algorithm that provides better performance in ver- 
tex reconstruction. We repeat the lifetime fit to this sample to determine the 
parameter values for R sig (At). We find the fraction / tai i for the tail part of 
the detector resolution is consistent with zero, and therefore set /tail = for 
this data sample 4 . In addition, the improvement in statistics enables us to 
determine some parameters that are previously determined only by MC sim- 
ulations. The fractions of the prompt component in R np for multiple- and 

• i i_ i i • / e multiple /-single r ~T)0 J ^multiple />singlc r r> — 

single-track vertices {J SB - , j s ^ tor ±r mesons, and J 5B - , J SB - for B 
mesons) are determined by the lifetime fit to the data. The values so obtained 
are consistent with the values determined using MC simulations. Table 6 lists 
the parameter values for i? sig (At). We find the resulting At resolution to be 
~ 1.43 ps (rms), improved over the resolution of ~ 1.56 ps obtained for the 
29 fb -1 sample. Figure 18 shows the observed At distributions for g£/ = +1 
(solid points) and g£y = —1 (open points) event samples. The asymmetry 
between the two distributions is proportional to sin 20 x and demonstrates the 
violation of CP symmetry. The value of sin 20! is measured to be 

sin 20i = 0.719 ± 0.074(stat) ± 0.035(syst). 

The quoted systematic error includes the systematic error due to the uncer- 
tainty in the resolution function (0.014). 



4 Consequently, s ta ii is not used 
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Table 6 

The resolution function parameters determined using a 78 fb _1 data sample. This 
data sample has been analyzed using a new track reconstruction algorithm. 



Parameters 


Values 


T B~o (P s ) 


1.551 ±0.018 


t b - (ps) 


1.658 ±0.016 



4,1 o.987i° : S 

sj ul 0.094 ± 0.008 



s 



t, 0.778 ± 0.048 



asc 



sL. 0.044 ± 0.002 



Smain 0.972 ± 0.045 

a ol (ps) 42.0lf;jj 

^multiple (1.65±S:g) X 10^ 

^single Q Q269 ± 00ig 

^multiple q ^ ± q Q42 



^multiple 44() ± 04 g 

/^ lc 0.701 ± 0.040 

ffz S - 0.764 ± 0.044 



6 Conclusions 



The resolution function, which is used in an unbinned maximum likelihood fit 
for the time-dependent measurements at the Belle experiment, is studied in 
detail. The resolution function is described as a convolution of three compo- 
nents; the detector resolution, the smearing due to non-primary tracks, and the 
kinematic approximation. The functional forms to describe these components 
are determined based on detailed MC simulation studies. The parameters for 
the detector resolution are determined using the data. The resulting resolution 
function has successfully described the At distribution used for measurements 
of B meson lifetimes, B°B° oscillation frequency, and the mixing- induced CP 
asymmetry parameter sin 20!. 
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